Universality in Heavy Fermions Revisited 
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A previous scaling analysis of pressure experiments in heavy fermion is reviewed and enlarged. 
We show that the critical exponents obtained from this analysis indicate that a one-parameter 
scaling describes these experiments. We obtain explicitly the enhancemente factors showing that 
these systems are indeed near criticality and that the scaling approach is appropriate. The physics 
responsible for the one-parameter scaling and breakdown of hyperscaling is clarified. We discuss a 
microsocopic theory that is in agreement with the experiments. The scaling theory is generalized 
for the case the shift and crossover exponents are different. The exponents governing the physical 
behavior along the non-Fermi liquid trajectory are obtained for this case. 
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I. INTRODUCTION 



The scaling theory of heavy fermions is based on the 
existence of a quantum critical point (QCP) which gov- 
erns the physical behavior of these systems |2|H • It yields 
the following scaling relations for the singular part of the 
free energy density 
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where the coherence temperature, T co h = \S\ UZ and the 
characteristic uniform field, h c = \5\^ h . H and h are 
staggered and uniform magnetic fields respectively. The 
scaling properties of relevant thermodynamic quantities, 
like the uniform susceptibility, xo oc d 2 f/dh 2 , the ther- 
mal mass, m T = C/T oc d 2 f/dT 2 , for T < T coh are 
obtained from Eq.l. The critical exponents obey stan- 
dard scaling relations but hyperscaling is modified due 
to the quantum character of the critical point. It is given 
by 2 — a = v(d + z) where z is the dynamic exponent 
and d the dimension of the system H. The quantity 
5 = (J/W) - (J/W) c measures the distance to the QCP. 
J and W are parameters of the Kondo lattice Hamil- 
tonian, the coupling between localized and conduction 
electrons and the bandwidth, respectively. 



II. SCALING ANALYSIS 

In heavy fermions the ratio (J/W) depends on volume 
V and consequently on pressure P. Let us define the crit- 
ical volume V c as the volume at which (J/W) = (J/W) c . 
Consider a physical quantity which close to V c behaves 
as X(P) = A\(V - V c )/V c \~ x where V is the volume at 
pressure P. If we introduce a reference pressure Po (vol- 
ume Vq) we have: 
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Denning, 5q = Vo — V c and AV = V — Vq, where the 
latter gives the change in volume due to the pressure 
change AP = P — Pq , we get 
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For AV sufficiently small, or small changes of pres- 
sure from the reference pressure, we have AV^ = 
—kqVqAP where the compressibility is given by kq = 
(-l/V)(dV/dP). In the limit that (AV/S ) < 1 we ob- 
tain 
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where (Vq/So) = Vo/(Vq — V c ) — ay /(ay — 1) with ay = 
Vo/V c . The equation above holds if the system at the 
reference pressure Po is not too close to the quantum 
critical point, otherwise the condition (AV/Sq) *C 1, or 
equivalently (kqAPolv) j '(1 — ay) <C 1 is not satisfied. 
Note that in Eq.5 the coefficient of AP depends on the 
critical exponent x associated with the physical quantity 
X. The validity of Eq. 5 for several physical quantities 
at and below the coherence line has been verified for the 
heavy fermions CeRu2Si2, CeAl^, UPt^ and CeCue as 
shown in Fig.l H. These materials are located in the 
non-critical side (V < V c ) of Doniach's phase diagram. 
For CeAl-j, a reference pressure of 1.2 kbars has been used 
to guarantee that this is the case, otherwise Po = 0. 
The collapse of the different data for a given material 
on a single line, as shown in Fig.l, implies the following 
relations among the critical exponents ||, 2 — a = vz 
and <ph — vz. The meaning and implications of these 
relations will be discussed further down. The inclination 



of the lines in Fig. 1, i.e., T\ 



Ln 



(2) different compounds are given in Table 1. 



/AP, for 
From Eq.5 
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we note that the Griineisen parameters fly = (Ty /kq) = 
xay /{ay — 1) provide essentially the enhancement factors 
due to the proximity of the quantum phase transition. 
Assuming, for example, that x = — 1 as we will discuss 
below, the results for fly m Table 1 yield values of ay 
ranging from ay ks 0.97 to ay w 0.99 as shown in this 
Table. This clearly indicates that the systems we are 
considering are close to the QCP and a scaling analysis 
is justified. 

Consider the case of CeCue. Eq.l yields, mj = C/T oc 
d 2 f/dT 2 oc \S\ ~ a ~ yz . From the experimental relation 
2-a = vz, we get mr' 1 oc \(V - V c )/V c \ vz = (l-ay)" 2 . 
Taking x = — 1, i.e., vz = 1, as before (see also below), 
we get an enhancement factor fly jay — 1/(1 — ay) for 
the thermal mass, m-r, of this material of approximately 
120. 

The data in Table 1 calls attention to the fact that 
Qy is larger for CeRu2Si2, although this system has the 
smallest thermal mass of the four compounds Q . In or- 
der to conciliate this result with the idea of universality, 
i.e., that the critical exponents are the same for a given 
quantity independent of the material, it is important to 
write the enhancement factor in terms of the parame- 
ters of the Kondo lattice Hamiltonian. For this purpose 
we need a relation between (J/W) and the volume. We 
assume that j(V) = (J/W) = (J/W) exp(-q(^^-)) 
where ( J/W)o is the value of this ratio at the reference 
pressure Pq or volume Vq and q a material-dependent 
parameter 0. For sufficiently small volume changes we 



get q(l — ay)/ay = (1 — aj) where aj 
(J/W) c = j(V c ). This allows us to write 
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where the parameter q, which relates changes in volume 
to changes in the interactions, depends on the particular 
system. Taking q — 5, for CeRu2Si2 and vz = 1, 
as before, we get an enhancement for the thermal mass 
of approximately 36 for this system. In this way taking 
different values of q for different systems we may explain 
the hierarchy of thermal masses using the same critical 
exponents. Note that for CeAl 3 the reference pressure 
Po = 1.2kbars such that mj- for this pressure is already 
very much reduced. 

Let us return now to discuss the scaling relations ob- 
tained from the data of Figure 1, namely, 2 — a = vz 
and <fih — vz. It is easy to verify that these equa- 
tions imply a simple one-parameter scaling such that 
Xq 1 oc m' T 1 oc h c oc A" 1 / 2 oc T co h, where A is the co- 
efficient of the T 2 term of the resistivity. As concerns 
the relations among the thermodynamic quantities they 
arise from a free energy obeying the simple scaling form, 
/ oc \5\ uz F[T /\5\ uz , h/\8\ vz }. This one-parameter scaling 
is reminiscent of single impurity and other phenomeno- 
logical, non-critical, approaches to the heavy fermion 
problem |jj . In the former case the characteristic temper- 
ature is identified with the Kondo temperature. There 



is however a fundamental difference between these ap- 
proaches and the present scaling theory. Here the char- 
acteristic energy scale given by the coherence tempera- 
ture T co h oc \{J/W) — (J/W) c \ uz vanishes at the criti- 
cal point of the Kondo lattice. Furthermore the scal- 
ing behavior found in these heavy fermion systems is 
due in our approach to their proximity to the QCP at 
(J/W) = (J/W) c , T = 0, H = 0, h = 0. The vanishing of 
at the QCP led to the prediction ||] and observation 
of non-Fermi liquid behavior in heavy fermion systems 

The observation of one-parameter scaling in a three di- 
mensional (d = 3) critical theory with three independent 
exponents is clearly associated here with the breakdown 
of the hyperscaling relation 2 — a = v(d + z). Note that 
the relation 2— a — vz arising from the experimental data 
is the hyperscaling relation for d = 0. It is not surprising 
then, that it yields scaling properties which are formally 
similar to those of a single impurity problem. Violation of 
hyperscaling is not uncommon in critical phenomena pj 
and below we shall discuss the possible reasons it occurs 
here. 

One-parameter scaling can also result from the con- 
straints imposed by the Fermi liquid (FL) behavior be- 
low the coherence temperature. Let us assume that the 
entropy of the local moments, of total angular momen- 
tum J, for T » Tcoh, S(T » Tcoh) = Nk B Ln(2J + 1) 
goes into that of the Fermi liquid which develops below 
Tcoh and is given by, S(T < T coh ) = f^° h dTC{T)/T 
with C(T) = ttitT '. Equating both entropies we find 
m^p 1 oc T co h, which in turn implies 2 — a = vz. This 
argument jl(^] relics on the existence of a single charac- 
teristic temperature in the non-critical side of the phase 
diagram which is not necessarily always the case, as will 
be seen below [0] . 



III. MICROSCOPIC MODEL 

We now discuss a model of nearly antiferromagnetic 
Fermi liquid |lj] which describes the results of the pres- 
sure experiments analyzed above. We start with the ex- 
pression for the singular part of the free energy density 
of a nearly antiferromagnetic Fermi liquid due to spin 
fluctuations O O lEl: 
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where 

X °(Q + q,u) = XQ(l-aq 2 -bLo 2 -Hew) (8) 
for q <C Q and u <C 1. At zero temperature we get 
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This can be rewritten as 
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where the correlation length £ = \5\ v and the distance 
to the QCP, (5 = 1 — I\Q- The correlation length expo- 
nent assumes the mean-field (or Gaussian) value, v — 1/2 
and the dynamic exponent, z = 2, typical of antifer- 
romagnetic spin fluctuations. All constants have been 
taken equal to 1. In the limit 6 — > 0, a change of 
variable yields the scaling properties of this model. We 
find / oc |<5| ly( - <i+z \ where d is the dimension of the sys- 
tem, with additional scaling corrections due to the term 
|5|(w£ z ) 2 . From Eq.l, with F[0,0,0] = constant we ob- 
tain the generalized hyperscaling relation, 2— a = v(d+z) 

§• 

Let us consider the case of localized spin fluctuations 
p2| , i.e., we neglect the q-dependence of x°(<Z,^)- We 
get 



f oc — / dui tan 1 
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Now for S — y 0, we find / oc \5\ vz , with scaling correc- 
tions. In this case, from Eq.l, we obtain 2 — a = vz 
and hyperscaling is violated due to the neglect of the q- 
dependence of x°. Note that in this local theory the ex- 
ponent vz = 1 as for the nearly antiferromagnetic Fermi 
liquid. 

For T ^ the relevant expression of the singular part 
of the free energy density is given by 

/ oc —T f dq ( dxcothxtan-^xTc^l- 



2 £2 



1 + q 2 e ' 

which can be written in the form 

/ oc \6\ v{d+z) F[T/T coh ] 



(12) 



(13) 



for x c = (3iv c /2 —1/2 and <? c £ very large, as close to the 
critical point. 

Again neglecting the q-dependent contribution we ob- 
tain, from Eq.12, the thermodynamic properties of the 
localized model. For xT^ z <C 1, i.e., w c £ z <C 1 we may 
write tan -1 y w y. Furthermore we take Tiu> c — ksT, in 
which case the constraint w c £ z <C 1 implies T <C T co h oc 
| S | .In this regime the free energy exhibits Fermi liq- 
uid behavior and is given by, / oc — (4/3)irq^g(l/2)T 2 £ i z , 
where g(y) = J dxx cothx. The breakdown of the scale 



invariant form, Eq.13 and consequently of hyperscaling 
is due to the neglect of the q-dependence. The thermal 
mass is given by 



-Wd^=m u M (l/3)<a*cx|«5| 



ic 1 (w) 



and diverges at the transition. m° M is a non-critical con- 
stant, which depends on local parameters and, linearly, 
on q c 

Within the same Fermi liquid regime, for T <C T co h, 
the free energy given by Eq.12 can be calculated and 
the thermal mass is given by, m^f = —d 2 f/dT 2 = 
m M [l — ((Jc() _1 ton _1 (g c ^)] . In this case the thermal 
mass increases but does not diverge as the system ap- 
proaches the critical point fljj. The second term in this 
expression for my is the universal scaling contribution, 
i.e., independent of the cut-off for q c ^ — > oc. It is pro- 
portional to oc l^ 1 ^ 2 since m° M oc q c . 

The free energies in the Fermi liquid regime (T <C T co h) 
as a function of the distance to the QCP, for the localized 
and nearly antiferromagnetic (NAF) models are shown 
in Fig. 2. Notice that for q c t; < 1/2, they nearly coin- 
cide. This can be seen directly, expanding my to obtain, 
Tot (<7c£ < 1) ~ m M (l/3)q 2 f = m£ (see Eq.14). Then, 
in this regime ]l6[ | the scaling properties of both models 
are the same and they yield similar results, for example, 
for the relative pressure variation of the thermal mass, 
m T (P )/m T (P). 

For the resistivity in the local approximation we find 
|7| 



Po 



dx- 
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For T£ z < 1, i.e., T < T coh = T z = \S\ V ', we get, 
p{T < T coh ) = P o(ir 2 /3){T/T coh ) 2 and for higher tem- 
peratures p(T » T coh ) = Po {n/2) {T/T coh ) |j. 

As concerns the magnetic field the experimental rela- 
tion (fih = vz implies that the characteristic field h c oc 
\S\ . In fact, in the local approach, the uniform mag- 
netic field h simply adds to the frequency a precession 
term Ml. Then at T = 



f oc - I dwtan" 1 \(w + h)£ z 

7T 



(16) 



and xo(T = 0) = d 2 f/dh 2 oc \5\- vz with vz = 1. This 
yields a diverging or enhanced uniform susceptibility but 
this result ceases to be valid sufficiently close to the crit- 
ical point (q c t; > 1). 

It is clear from the results above that the theory 
of antiferromagnetic local spin fluctuations, with a sin- 
gle characteristic energy scale, T coh oc \S\ UZ , such that 
f(T = 0) oc \Sr = T coh , m T oc \8\ 



-VZ rp— 1 



T~ oh and p(T < T coh ) 
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Xo (T = 0) oc \5\-» z 
with A oc T~ 2 h correctly describes the pressure experi- 
ments we analyzed before. Furthermore in this approach 
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the crossover exponent vz = 1. This provided the moti- 
vation for assuming this value for vz in the calculation 
of the enhancement factors in Section II of this paper. 

The theory of localized antiferromagnetic param- 
agnons can be summarized in the scaling form of the 
free energy, / oc \5\ VZ F[T/\5\ VZ , h/\5\ vz ] with vz = 1. Al- 
though this is in agreement with the pressure experiments 
in a region of the phase diagram, for ( J/W) > ( J/W) c , 
T < T co h, as the systems get closer to the critical point 
and q c (; — > oo the full g-dependence of the dynamic sus- 
ceptibility must be taken into account. The local the- 
ory also does not describe the non-Fermi liquid behavior 
observed at {J/W) = {J/W) c , T -> 0. For example, 
it's result for the specific heat at |<5| — is given by, 
C = (3/2)NkB, i.e., that of a classical gas of N free par- 
ticles. Sufficiently close to the QCP it is then necessary 
to consider the g-dependence of \° (Q + 1, w ) • Taking this 
into account we obtain the universal contribution for the 
free energy at the QCP, f[S = 0] oc T~z~ , in agree- 
ment with Eq.13, and the specific heat C/T(8 = 0) oc 
d 2 f/dT 2 oc T 1 / 2 for d = 3, z = 2 with vz = 1 PJMH. 



and one finds the previous result, C/T oc . However 
in the situation here, where ip = 2/3 < vz = 1, the tilde 
exponents also play a role in determining the behavior 
along this path, which is tangent to the Neel line at the 
quantum multicritical point as shown in Fig. 3. We find 
for the specific heat 



C/T oc T~ 



(18) 



Assuming thermal Gaussian exponents, essentially a — 
1/2, wc get, C/T oc T 5 / 4 for tp = 2/3, v = 1/2 and z = 2, 
instead of C/T oc T 1 / 2 for the case of extended scaling. 
The staggered susceptibility xq{5 = 0,T) oc T"^ = 
T-z/ 2 since 7 = 7= 1 @. 

An interesting possibility has been raised by Rosch 
at al. Q who claimed that two-dimensional fluctua- 
tions are those relevant to describe the observed crit- 
ical behavior. In this case scaling is extended since 
if) = z/(d + z — 2) = vz = 1. This approach leads to 
C/T(S = 0) oc LnT at the QCP and in the Fermi liquid 
regime, for T <C T co h, to tut oc Ln\S\. 



IV. GENERALIZED SCALING 

It is important to point out that the result C/T oc T 1 / 2 
at I (5 1 — which arises from the scaling form Eq.13 is valid 
only in the case of extended scaling, i.e., if) = vz. Here 
if) is the shift exponent such that the Neel temperature, 
T N oc |<5|^, close to the QCP MO (see Fig.3). It turns out, 
experimentally, that ip = 1 |9fhowever, the oretically , one 
z/(d + z 



19,11. In 



gets, ip = z/{d + z - 2) = 2/3 ^ vz = 1 |4Jj_ 
this case a more general scaling form for the free energy 
is required to describe the complex critical behavior in 
the neighborhood of the QCP Q. It is given by 0], 

/oc \S(T)\ 2 - a F c [t] 



T 



\S(T)\ 



(17) 



with 



5(T) = S(T = 0) - T 1 ^ 



The singularities along the Neel line, |<$(T)| = 0, are 
described by a new set of tilde exponents a, v, etc., 
different from those associated with the zero tempera- 
ture fixed point (the non-tilde exponents). The scaling 
function F c [t = 0] = constant and F c [t — > 00] oc t x 
with x — (a — a)/vz such that close to the critical 
Neel line we obtain the correct asymptotic behavior, 
/ oc A(T)\5(T)\ 2 ~ a , where the amplitude A(T) = T^r . 
It is easy to verify that in the case of extended scaling, 
i.e., ip = vz, the exponents associated with the zero tem- 
perature fixed point are sufficient to describe the behav- 
ior along the non- Fermi liquid trajectory, \5\ = Q, T — *■ 



V. ABOVE THE UPPER CRITICAL DIMENSION 

In quantum phase transitions the relevant dimension- 
ality is d e ff — d + z, as is evident from the modified hy- 
perscaling relation, 2 — a — v(d + z). In the problems we 
have studied above, it turns out that d e ff > d c , where, 
d c = 4 is the upper critical dimension for these magnetic 
transitions. This is the reason Gaussian theories, as the 
SCR theory of spin fluctuations Jl^ | , provide an adequate 
description of the quantum critical point in nearly ferro 
and antiferromagnetic 3d materials ]l9| . Let us consider 
further implications of the fact that d + z > d c = 4. Con- 
sider the expression for the singular part of the T = 



free energy, / oc \S\ 



v(d+z) 



Since the Gaussian expo- 



nent v = 1/2, whenever d + z > 4, we can rewrite it 
as, / oc \5\ 01 with a < 0. In this case an analytic ex- 
pansion of the free energy close to the critical point, such 
that, / oc \S\ 2 , will always dominate the Gaussian con- 
tribution for S sufficiently small. The total free energy, 
in the non-critical side of the phase diagram, below the 
coherence line, can be written as a sum of an analytic 
term and a Gaussian contribution, 



ft = \6\ 2 + \5\ v{d+z) F\ 



T 

w 



(19) 



where we neglected the temperature dependence of the 
analytic part assuming it is less singular than that of the 
Gaussian. In the Fermi liquid regime, i.e., T -c T co h — 
\5\ vz , we have ||, 



ft = \Sf + \S\ 



v(d+z) 



T 

w 



(20) 



then, sufficiently close to the critical point and for d+z > 
4, with v = 1/2, we obtain, 



4 



ft = \5\ 2 + \5\» (d+z) (^) 2 + -" ( 21 ) 
which can be rewritten in the scaling form 

/* = \5\ 2 F[^-\ (22) 

where the value of a = in the equation above, see Eq.l, 
can be associated with the breakdown of hyperscaling for 
d+z > 4. The new spin-fluctuation temperature is given 
by 

Tsf = (23) 

where we used v = 1/2. Consequently the effect of the 
analytic contribution is to introduce a new energy scale in 
the Fermi liquid region of the phase diagram [ plj , namely 
T s f. If we calculate the specific heat from the above 
expression for the free energy, we get, C/T oc |(5|^ d z ' 
which of course coincides with the Gaussian result. On 
the other hand taking into account the field depen- 
dence of the analytic (mean-field) and Gaussian contri- 
butions, it can be easily shown that the order param- 
eter linear susceptibility, xo = (d 2 f j 'dh 2 ) with / oc 

\5\ 2 F (T/T sf ,h/\6f + ~<) is given by xo = \5\- l F x {T /T sf ), 
for T <C T co h 1 where we used the mean field exponents, 
(3 = 1/2 and 7=1. The field h here is that conjugated 
to the order parameter. 

We point out also that, because d + z > 4, even at 
\S\ = there is a characteristic field h cross — J which 
yields the crossover from mean-field to Gaussian behavior 
p| . J is the actual critical coupling between localized 
and itinerant electrons Q . The order parameter m at the 
quantum critical point varies with the conjugated field 
h, according to m oc (h/J) 1 / 5 [Q. For small fields, i.e., 
h <C h cross , the mean- field contribution with Smf = 3 
dominates while in the opposite case the Gaussian one 
with Sg = 01 ^ s dominant (in the situation of 

interest here d = 3 and z — 2 such that Sg = 7/3 < 
Smf = 3). 

Finally note that the above discussion does not affect 
the local spin fluctuation results since, in this case the 
T = free energy is more singular than the analytic 
contribution, at least for vz < 2. 

VI. CONCLUSIONS 

Our analysis of the pressure dependence of several 
physical quantities for different heavy fermions on a re- 
gion of the phase diagram, for {J/W) > (J/W) c , T < 
T^, has shown that these systems are close to a quan- 
tum critical point and a scaling approach is indeed ap- 
propriate. We have obtained from the experimental data, 
Griineisen parameters which yield the enhancement fac- 
tors due to the proximity of the QCP. It turns out that 



a one-parameter scaling is sufficient to describe the ex- 
periments in this region of the phase diagram and this is 
associated with the violation of the hyperscaling relation. 
This one-parameter scaling here is different from that of 
single impurity approaches where the characteristic en- 
ergy scale is identified with the Kondo temperature. In 
our case T co h vanishes at the quantum critical point and 
this leads to the appearance of non- Fermi liquid behavior 
at S = 0. We have shown that a theory of localized spin 
fluctuations, with vz = 1, describes the pressure exper- 
iments summarized in Fig.l. For these experiments the 
q-dependence of the dynamic susceptibility plays no role. 
Physically this must be related to the fact that the actual 
spectrum of spin fluctuations is nearly ^-independent, at 
least in the relevant directions of g-space, in this part of 
the phase diagram. This region corresponds to the case 
<7c£ < 1/2 of the NAF model since, as we have shown, 
the scaling properties of this and the local model are the 
same in this regime. In fact the local approach can be ob- 
tained as an expansion of the full g-dependent theory of 
antiferromagnetic spin fluctuations for q c £ <C 1 (see also 
Jl6[). Then it holds in a region, not too close to the QCP, 
but where scaling, specific to this q-independent regime 
(one-parameter scaling), still applies. The local theory 
allows for an explicit calculation of the Wilson (rar/xo) 
and Kadowaki- Woods {A/m\) ratios Jj| which turn out 
to be constants, i.e., independent of the distance to the 
critical point. 

As the system moves closer to the critical point (q c £ — » 
oo) we leave the limit of validity of the local theory and 
the full q-dependent susceptibility must be used. In both 
cases vz = 1. In the q-dependent regime a generalized 
scaling theory is required to describe the complex critical 
behavior in the neighborhood of the QCP since the shift 
and crossover exponents are different for nearly antiferro- 
magnetic 3c? systems. This is particularly relevant along 
the non-Fermi liquid trajectory which is tangent to the 
QCP for t/j = 2/3 < 1. It turns out that the exponents 
along this line depend on the thermal exponents charac- 
terizing the singularities on the Neel line, besides those 
associated with the QCP. 

Although for the systems investigated the effective di- 
mension d + z > d c , it is the local nature of the spin 
fluctuation spectrum in the part of the phase diagram 
investigated that is responsible for the breakdown of hy- 
perscaling and one-parameter scaling, observed in the ex- 
periments. In the microscopic approach the violation of 
hyperscaling can be traced to the cut-off acting as a rel- 
evant variable. 

The purpose of a scaling theory of heavy fermions is 
to provide a unified description, in terms of a set of criti- 
cal exponents, not only of the non-Fermi liquid regime 
but also of both sides of the quantum critical point. 
While most of the studies have been carried on in the 
non-critical side we also expect to find scaling behav- 
ior on the ordered region of the phase diagram, i.e., for 
( J/W) < (J/W) c , where the systems at T = have long 
range magnetic order. We hope future studies will also 
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include this interesting region and also a larger class of 
materials [E2fl. 
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APPENDIX: 



AT 2 f 2 

%(1/% C 3 



J Q a? 



12AT 2 tf , , . , 



0(A 2 ) 



(28) 



The hrst term is independent of A and corresponds to 
the localized model where the stiffness A vanishes. The 
validity of the local theory is then given by q c £, < 1 or 
? c £l *C (Jqo. 2 /A) 1 / 2 such that when A — ► the local 
model becomes valid arbitrarily close to the QCP. 



We discuss here the limit q c l; <§; 1 in relation to the 
theory of Ref. @. The free energy is given by, 



/=- 



3T 



dq 



dx coth x tan 



xT 



Aq 2 
(24) 



where x = Jq is the q-dependent exchange coupling 
between the local moments at the wavevector Q of the in- 
cipient magnetic instability and Jq its critical value. A is 
the stiffness of the spin fluctuations defined by the small 
wavevector expansion of the magnetic coupling close to 

the wavevector Q, i.e., Jq — JQ+ q = Aq 2 H . 

A redefinition of the relevant quantities yields, 



s = J t 



dq 



dx coth x tan 



xTC 



A(l + q 2 e) 



(25) 



where the correlation length £ may be written as, £ = a^L 
with a = {A/ Jqo 2 ) 1 / 2 , & = a^" 1 / 2 and 5 = 1- 
( Jq/Jq), such that, v = 1/2. a is an interatomic distance 
and the dynamic exponent z = 2. Since 1/(1 + y 2 ) < 1, 
for XcT ^- <C 1, or equivalently, w c ££ <C 1, which assum- 
ing TiLo c = ksT implies T <C r co ^ oc £7 Z , we can expand 
the tan -1 for small arguments. In this regime the free 
energy exhibits Fermi liquid behavior. Performing the 
integrations, we get 



12T 2 ^- rf ) / tan- 1 ^ 

/ ~ \ ff(l/2)?cC 1 - 



.4 



where g(y) = x cothxdx. The limit q c £ <gc 1 is ob- 
tained, either because the system is not too close to the 
QCP, i.e., \S\ is large, or A — > due to the localization 
of the fluctuations, we get 



12T 2 
/ cx —g(l/2)q c 



Using £ — a^L, as defined above, we get 
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TABLE 



Compound 


Tv 


K 


Q.v — Fv / Ko 


ctv = Vq/Vc 


CeAh 


89 


2.17 


41 


0.976 


CeRu2Si2 


171 


0.95 


180 


0.994 


UPt 3 


26 


0.48 


54 


0.981 


CeCu 6 


133 


1.1 


121 


0.992 



TABLE I. Griineisen parameters for different heavy 
fermions according to the relative pressure variation of several 
physical quantities shown in Fig. 1 Q . Ty and the compress- 
ibility kq are in Mbar~ x . The reference pressure for CeAlz 
is Po = 1-2 kbars otherwise Po = 0. The data for C'eCue is 
taken from Ref. [HI and references therein. 



FIGURE CAPTIONS 

Figure 1. Semi-logarithmic plot of X(P)/X(P ) for 
several physical quantities, at or below T co h, as a function 
of pressure for different heavy fermions |3|. For CeAl^, 
Pa = 1.2 kbars otherwise Po = 0. The numbers close to 
the lines are their inclinations, Ty, given in Table I (also 
see text). 



Figure 2. The free energy in the Fermi liquid regime, 
at a fixed temperature To <C T co h, for the localized and 
nearly antiferromagnetic models as a function of the dis- 
tance to the critical point, located at (<7c£) -1 = 0. For 
(icO^ 1 > ^ the free energies of both models nearly coin- 
cide and consequently have the same scaling behavior. 



Figure 3. Phase diagram of a nearly antiferromagnetic 
3d system, with the shift exponent if) — 2/3 and the 
crossover exponent vz = 1. The tilde (thermal) expo- 
nents determine the critical behavior on the Neel line. 
In this case of ip < 1, the non- Fermi liquid trajectory is 
tangent to the quantum critical point and the thermal 
exponents besides those associated with the QCP are re- 
quired to characterize the physical behavior along this 
line (\5\ = 0,T->0). 
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Pressure (kbars) 
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